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Popular Summary 

During the months of July- August 2007 NASA conducted a research campaign called the 
Tropical Composition, Clouds and Climate Coupling (TC4) experiment. Vertical profiles of 
ozone were measured daily using an instrument known as an ozonesonde, which is attached to a 
weather balloon and launch to altitudes in excess of 30 km. These ozone profiles were measured 
over coastal Las Tablas, Panama (7.8N, 80W) and several times per week at Alajuela, Costa Rica 
(10N, 84W). Meteorological systems in the form of waves, detected most prominently in 100- 
300 m thick ozone layer in the tropical tropopause layer, occurred in 50% (Las Tablas) and 40% 
(Alajuela) of the soundings. These layers, associated with vertical displacements and classified 
as gravity waves (“GW,” possibly Kelvin waves), occur with similar structure and frequency 
over the Paramaribo (5.8N, 55W) and San Cristobal (0.92S, 90W) sites of the Southern 
Hemisphere Additional Ozonesondes (SHADOZ) network. The gravity wave labeled layers in 
individual soundings correspond to cloud outflow as indicated by the tracers measured from the 
NASA DC-8 and other aircraft data, confirming convective initiation of equatorial waves. 

Layers representing quasi-horizontal displacements, referred to as Rossby waves, are robust 
features in soundings from 23 July to 5 August. The features associated with Rossby waves 
correspond to extra-tropical influence, possibly stratospheric, and sometimes to pollution 
transport. Comparison of Las Tablas and Alajuela ozone budgets with 1999-2007 Paramaribo 
and San Cristobal soundings shows that TC4 is typical of climatology for the equatorial 
Americas. Overall during TC4, convection and associated meteorological waves appear to 
dominate ozone transport in the tropical tropopause layer. 
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35 Convective and wave signatures in ozone profiles over the equatorial 

36 Americas: Views from TC4 (2007) and SHADOZ 

37 Abstract. During the TC4 (Tropical Composition, Clouds and Climate Coupling) 

38 campaign in July-August 2007, daily ozonesondes were launched over coastal Las Tablas, 

39 Panama (7.8N, 80W) and several times per week at Alajuela, Costa Rica (10N, 84W). 

40 Wave activity, detected most prominently in 100-300 m thick ozone laminae in the 

41 tropical tropopause layer, occurred in 50% (Las Tablas) and 40% (Alajuela) of the 

42 soundings. These layers, associated with vertical displacements and classified as gravity 

43 waves (“GW,” possibly Kelvin waves) by laminar identification, occur with similar 

44 structure and frequency over the Paramaribo (5.8N, 55W) and San Cristobal (0.92S, 90W) 

45 Southern Hemisphere Additional Ozonesondes (SHADOZ) sites. GW-labeled laminae in 

46 individual soundings correspond to cloud outflow as indicated by DC-8 tracers and other 

47 aircraft data, confirming convective initiation of equatorial waves. Layers representing 

48 quasi-horizontal displacements, referred to as Rossby waves by the laminar technique, are 

49 robust features in soundings from 23 July to 5 August. The features associated with 

50 Rossby waves correspond to extra-tropical influence, possibly stratospheric, and 

5 1 sometimes to pollution transport. Comparison of Las Tablas and Alajuela ozone budgets 

52 with 1999-2007 Paramaribo and San Cristobal soundings shows that TC4 is typical of 

53 climatology for the equatorial Americas. Overall during TC4, convection and associated 

54 waves appear to dominate ozone transport in the tropical tropopause layer; intrusions 

55 from the extra-tropics occur throughout the free troposphere. 
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56 l. Introduction 

57 Ozone in the tropical troposphere reflects an interaction of photochemical and 

58 dynamical factors. The marine atmosphere is usually unpolluted, largely because the 

59 boundary layer (BL; see Appendix for abbreviations) is a region of 0 3 photochemical loss 

60 [ Piotrowicz et al., 1991]. This is a consequence of slow formation (low NO x conditions; 

61 Zafiriou et al., 1980; Thompson et al., 1993) with loss by HO x , or occasionally by rapid 0 3 

62 destruction by halogens [ Read et al., 2008]. In the free troposphere mixed sources 

63 converge [Thompson et al., 1996]. Advected pollution, stratospherically influenced air, 

64 and lightning add to 0 3 formation, the latter at rates according to time since a lightning 

65 episode ( Thompson et al, 1997; see Cooper et al. [2006] and Bertram et al. [2007] for 

66 analyses of lightning influence in mid-latitudes). Extra-tropical ozone may also enrich 

67 ozone in the tropics [ Randel et al., 2007]. 

68 Examination of ozone profiles from sondes or aircraft over remote tropical sites 

69 reveals the free troposphere as a region of low 0 3 (< 30 ppbv; Thompson et al., 2003a) 

70 alternating with layers of elevated 0 3 (sometimes > 100 ppbv; Newell et al., 1999). From 

71 soundings taken through SHADOZ (Southern Hemisphere Additional Ozonesondes; 

72 Thompson et al., 2003a, b; 2010; submitted article available as Supplemental Material), 

73 SOWER [Stratospheric Ozone and Water in Equatorial Regions; Hasebe et al., 2007; 

74 Takashima and Shiotani, 2007] and related campaigns, much is known about the 

75 structure of 0 3 in the tropical upper troposphere (UT) and lower stratosphere (LS). Stable 

76 pollution layers in the free troposphere are frequent over Reunion, Fiji, Samoa, San 

77 Cristobal, and Ascension [ Thompson et al., 1996; 2003b; Oltmans et al., 2001; 2004; 

78 Randriambelo et al., 2003]. Reduced O., layers often characterize the UT and TTL 

79 (tropical tropopause layer, -8-14 km; Fueglistaler et al., 2009), where convective outflow 

80 of low-0 3 BL air occurs. If 0 3 in the UT or TTL averages to lower concentrations than in 

81 the mid-troposphere, an “S-shape” profile results [ Folkins et al., 2000]. 

82 The laminar identification (LID) method, based on 0 3 and potential temperature 

83 gradients [Teitelbaum et al., 1994, 1996; Grant et al, 1998; Thompson et al., 2007a; 

84 2008], interprets persistent 0 3 layers in terms of two wave-types. In the tropics, Rossby 

85 waves (RW) represent horizontal displacement; these correspond to filaments of extra- 

86 tropical air or pollution from long-range transport. When SHADOZ profiles are analyzed 

87 with the LID technique, RW signatures appear in < 20% of the soundings [ Loucks , 2007; 

88 Thompson et al., 2010]. Signatures of convectively-generated gravity waves (GW) occur 
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89 in 40-90% of SHADOZ soundings, depending on location and season [ Thompson et al., 

90 2010a, b]. Near the tropopause, GW are usually identified with Kelvin waves. Transient 

91 Kelvin waves have been observed in 0 3 over both eastern and western Pacific [ Fujiwara et 

92 al., 1998; 2001]. 

93 Wave activity over the equatorial Americas has received less attention. Robust GW 

94 and RW signals were noted during the Milagro/INTEX-B/IONS-06 (Intercontinental 

95 Transport Experiment; INTEX Ozonesonde Network Study) campaigns over Mexico City 

96 (19N, 99W) and Houston (30N, 95W) in March-May 2006 [ Thompson et al., 2008], two 

97 locations that are essentially sub-tropical when springtime air parcel flows link them to 

98 central America [ Fast et al., 2007]. The TC4 (Tropical Composition, Cloud, and Climate 

99 Coupling) mission in July- August 2007 offered an opportunity to characterize 0 3 profiles 

100 closer to the equator than the IONS-06 soundings. TC4 [ Toon et al., 2010] investigated 

101 chemical transformation in convective systems and the impacts of deep convection on 

102 constituent transport, dehydration and cirrus formation. Sampling from San Jose, Costa 

103 Rica, with NASA’s DC-8, WB-57 and ER-2 aircraft, was well-suited for comparisons with 

104 observations from Las Tablas, Panama (LTP, 7.8N, 80W), where daily soundings were 

105 made from the NATIVE (Nittany Atmospheric Trailer and Integrated Validation 

106 Experiment) mobile lab. Most TC4 flights were south of the Intertropical Convergence 

107 Zone (ITCZ), which was located at 12-13N during the experiment [ Toon et al., 2010]. 

108 SHADOZ [Thompson et al., 2003a] soundings from Costa Rica (Alajuela, ioN, 84W) were 

109 augmented with several launches per week from early July through 9 August 2007. 

1 10 Five TC4 studies are devoted to analysis of ozone transport and wave activity. In 

1 1 1 Selkirk et al. [2010] temperature anomalies from four-times-daily radiosondes over 

1 12 Alajuela indicate equatorial waves in the TTL. Comparison of TC4 radiosondes and water 

1 13 vapor from cryogenic frost-point hygrometer readings suggests less convection than 

1 14 during the 2005 TCSP (Tropical Convective Systems and Processes) campaign in the 

1 15 same location and time of year. The free troposphere was drier in TC4 than TCSP, and 0 3 

1 16 in the mid-upper troposphere was ~35% higher in TC4 (Figures 1, 3 and 4 in Selkirk et 

1 17 al., 2010). Avery et al. [2010] focus on the strong anti-correlation between ozone and 

1 18 condensed cloud water content measured during sampling in active convection, and 

1 19 generalize their findings on convective redistribution of ozone with tracer composites. 

120 Convective outflow as indicated by local 0 3 minima and elevated lower tropospheric 

121 tracers appears to maximize at 10-11 km. Avery et al. [2010] argue that in the DC-8 


4 



122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 


sampling region near the ITCZ, there has been about 50 % convective turnover below the 
TTL, with vertical transport from just above the boundary layer. In Morris et al. [ 2010 ] 
lower-mid tropospheric 0 3 within a convectively active region over the Pacific southwest 
of the Las Tablas (LTP) launch site on 5 August 2007 is investigated in detail. The LTP 
sounding for that day was caught in a strongly convective cell that bounced the balloon 
package up and down between 2 . 5-5 km; the 0 3 within this layer increased from 30 to 40 
ppbv before the balloon ascended to the stratosphere. Petropavlovskikh et al. [ 2010 ], 
analyzing the 17 July 2007 DC -8 flight, find that influences on 0 3 near the TTL are a 
mixture of convection and advection, the latter including extra-tropical air, some of it 
possibly of stratospheric origins. The lowest 0 3 layer occurred at 9-11 km, but low-0 3 
signatures were also recorded by the UV-DIAL instrument at 2-4 km and 13 km. 

As insightful as the above investigations are, none systematically examines 
dynamic influences in ozone structure in the free troposphere and LS day-by-day during 
TC 4 . That is the purpose of the present study. Using the full record of Alajuela (ACR) 
and Las Tablas 0 3 soundings (Section 2), the following analyses are performed: 

1 ) Mean 0 3 profiles are determined and day-to-day variability is examined in 
curtains of 0 3 mixing ratio and in integrated 0 3 column amounts. 

2 ) Ozone budgets based on LID and expressed as amounts affected by GW and 
RW provide a consistent framework for examining dynamic influences over the course of 
TC 4 . The climatology of wave signatures at LTP and ACR is compared to those at 
Paramaribo ( 5 - 8 N, 55 W) and San Cristobal (o. 92 S, 9 oW) sondes, long-term SHADOZ 
stations also in the equatorial Americas (Section 3.1). 

3 ) Case studies of soundings and ancillary aircraft measurements are used to 
corroborate wave designations, with convection for GW and stratospheric influences or 
pollution for RW (Section 3.2). One episode is among those analyzed by Avery et al. 
[ 2010 ] and Morris et al. [ 2010 ]. Others were chosen to illustrate contrasting impacts. 

4 ) In Section 4 context for the TC 4 observations is given by June-July-August 
(JJA) Costa Rican sondes in 2006 , and the 9 -year SHADOZ record [Thompson et al., 
2003 a, b] at Paramaribo and San Cristobal. Specifically we ask: 

> How does tropospheric 0 3 over LTP and ACR in 2007 compare to JJA 0 3 over 
Costa Rica (Heredia site near Alajuela) in 2006 , when an El Nino affected 
the eastern Pacific [Arguez, 2007 ]? How do LTP and ACR 0 3 during TC 4 
compare to 0 3 over San Cristobal in 2007 ? 
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> How does 0 3 over the equatorial Americas in TC4 compare to the 1999-2006 
record? Indices based on GW activity are used to quantify interannual 
variability. 

2. Experimental. Observations and Methods of Analysis. 

2.1 Ground & Aircraft 

Continuous surface 0 3 measurements at Las Tablas (7.8N,8oW) were made during 
the period 13 July to 10 August 2007 from the NATIVE mobile lab that also included an 
ozonesonde ground station. The surface 0, } was measured with a TECO Model 49 C ozone 
analyzer, along with carbon monoxide (TECO Model 48CTL), NO and NO y (TECO Model 
42CY), S0 2 (TECO Model 43C-TLE) and particle size distribution (Scanning Mobility 
Particle Sizer). All data can be viewed at < http://ozone.met.psu.edu/NATIVE/ 
measurements. html>. Calibrations of 0 3 , CO, and S0 2 were made prior to and directly 
after the campaign with instrument grade gases (Airgas, Inc.). Calibration of NO and NO y 
was performed daily with instrument grade NO (Airgas, Inc.). Catalytic conversion 
efficiency tested before and after TC4 remained close to 100%. 

2.2 Ozone Profiles and P-T-U Profiles 

All ozone profile data analyzed here were taken with electrochemical concentration 
cell (ECC) instruments coupled with standard radiosondes, as described in Thompson et 
al. [2003; 2007a]. Coordinates of LTP, ACR and the SHADOZ sites referred to here 
appear in Table 1, along with details of the sondes used. At Las Tablas, the 0.5% KI 
buffered solution with ENSCI ozonesondes is a combination that optimizes tropospheric 
and stratospheric 0 3 measurements [Smit et al., 2007; Thompson et al., 2007b; Deshler 
et al., 2008]. Vaisala radiosondes, Model RS-80, were used to collect P-T-U data. For 
nine of the 25 LTP soundings during TC4 the humidity data are unreliable above ~300- 
500 hPa due to suspected sensor icing. Ozone profiles from LTP are viewable at 
< http://ozone. met.psu.edu/Panama Data/index.html>. Each ozonesonde is calibrated 
according to standard procedures prior to launch; all LTP sondes were also compared to 
the TECO 0 3 for 5-10 minutes pre-launch; agreement is within the stated precision of 
each technique (5%; Figure 1 in Morris et al., 2010). 

Ozone over Alajuela was measured by ENSCI ECC sondes with 1% KI with reduced 
(0.1%) buffer (Table 1); RS-80 radiosondes with a cryogenic frost-point hygrometer were 
used for P-T-U. At both ACR and LTP, most launches took place in early afternoon local 
time to capture Aura [ Schoeberl et al., 2006], Aqua and CALIPSO satellite overpasses 
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[Toon et al., 2010]. Images and data for all Costa Rican (late 2005 start), Paramaribo 
(1999-present) and San Cristobal (1999-present) ozonesonde and P-T-U profiles are 
available at < http://croc.gsfc.nasa. gov/shadoz> and at the World Ozone and Ultraviolet 
Data Centre, < http:/ /woude.org> . 

Aircraft data used most often in analysis of the soundings are: 0 3 from the 
FASTOZ in-situ instrument [ Avery et al., 2010], uv-DIAL and DACOM CO on the DC-8; 
the Cloud Physics Lidar (CPL) and Cloud Radar System (CRS) on the ER-2 [ McGill et al., 
2004; Hlavka et al., 2010]. Regional cloud and convective information comes from 
meteorological analyses and GOES imagery, as archived by Toon et al. [2010]. 

2.3 Ancillary Data 

As for IONS-04 [ Thompson et al., 2007a] and IONS-06 [Thompson et al., 2008], 
tracers for 0 3 origins include: (1) RH from the radiosonde P-T-U profiles; (2) Ertel’s 
potential vorticity (pv; 1 pvu = to' 6 m 2 s-i/K) computed from the Goddard Earth 
Observing System Assimilation Model (GEOS-version 4; Bloom et al., 2005); (3) forward 
and backward air-parcel trajectories for each launch location and date, calculated with 
the kinematic version of the GSFC trajectory model [Schoeberl and Sparling, 1995] using 
GEOS meteorological fields at a lxi-degree grid. Lightning imagery is also used to 
describe potential 0 3 influences, as are absorbing aerosol data from OMI, trajectory- 
enhanced aerosol-exposure images and OMI N 0 2 amounts. All back-trajectories in the 
TC4 region, aerosol and lightning data and trajectory-mapped exposure products are at 
the GSFC TC4 website: < http://croc.gsfc.nasa.gov/tc4>- In selected cases, additional 
trajectories were run with the NOAA Hysplit model ( Draxler and Rolph, 2003). 

2.4 Analysis for Wave Influences 

In the LID (Laminar Identification; Thompson et al., 2007a; 2008) technique, 0 3 
and potential temperature ( 0 ) laminae, as described in Teitelbaum et al. [1994; 1996] and 
Pierce and Grant [1998] are used to identify RW and GW signatures as illustrated in 
Figure 1. The method is described as follows: 

1) For each sounding, the 0 3 and 0 laminae are isolated through normalization to 

running means. A boxcar smoothing over 2.5 km is used to isolate each lamina. 
Larger and smaller ranges (from 0.5 to 2.5 km) have been tested with little 
difference in laminar statistics. Normalized 0 3 is the solid line in Figure 1; the 0 
deviations are signified by the dotted line. Because the ozone mixing ratio 
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precision is 5 % in most of the region of interest, only laminae representing a 
deviation of _> 10 % are included in the analysis (red vertical lines in Figure 1). 

2 ): Correlations between 0 3 and 0 gradients are compared to identify GW and RW. 

Where the correlations (as in dashed line, Figure 1) exceed 0.7 GW (light green) 
is defined, based on the reasoning that such correlation between the two quantities 
signifies a vertical disturbance. Conversely, when horizontal motion affects 0 3 , it 
is poorly correlated with 0 and RW is specified. We adopt the practice of Pierce 
and Grant [ 1998 ] and Teitelbaum et al., 1994 ], setting the RW limit for anti- 
correlation at + 0 . 3 , corresponding to light blue sections in Figure 1. This is a 
fairly conservative criterion; see Thompson et al. [ 2007 a] for sensitivity studies. 

Two analyses are performed with the LID results. First, for an ensemble of soundings, 
wave frequencies at a given altitude are calculated from the percentage within a sample set 
that have laminae designated as RW or GW. Second, for each sounding, the contributions 
of RW and GW above the boundary layer (BL) to the tropopause or to 20 km are 
computed by integrating the amount of 0 3 within a given layer and adding up all the RW 
and GW segments. Thus, “tropospheric ozone” amounts in this paper are actually free 
tropospheric column amounts, above a BL top that varies from 1.2 km (San Cristobal) to 
1.8 km (Paramaribo). The BL top is determined by taking the most negative second 
derivative of 0 with respect to altitude, between 0.4 and 2.5 km [Yorks et al., 2009 ]. 
Dobson Units (DU) are used; one DU = 2.69 x 10 16 cm' 2 . The amount of 0 3 within the 
column not identified with RW or GW is labeled “other” in the budgets. 

For determination of tropospheric LID 0 3 budgets (Section 3 ), an ozonopause 
definition of tropopause is employed (white dots, Figure 2; using the method of Browell 
et al., 1996 ). This technique, in which 0 3 gradients are approached from the stratosphere, 
falls closer to the cold-point tropopause than the ozonopause of Selkirk et al. ( 2010 ; their 
Figure 2 , from Alajuela soundings). Thompson et al. [ 2007 a] showed that tropospheric 0 3 
columns can differ significantly under certain conditions, depending on the use of 0 3 or 
thermal tropopause. Such occurrences are infrequent in mid-latitudes (< 10 % in IONS -04 
or IONS -06 soundings) and are negligible in TC 4 , where Selkirk et al. ( 2010 ) showed that 
the thermal tropopause and chemical tropopause, based on either water vapor from frost- 
point hygrometry or 0 3 mixing ratios, was consistently at 350-355 K (147 hPa, 14.2 km); 
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the cold-point tropopause was always 370-378K (16-17 km). For a 1-2 km difference in 
thermal tropopause and the higher ozonopause, the 0 3 budget might vary 2-3 DU. 

3 . Results and Discussion: July-August 2007 
3.1. Overview of Ozone Profiles over LTP and ACR 

From time-series of 0 3 mixing ratio over Las Tablas and Alajuela (Figure 2) the 
following features emerge: (1) Roughly 10% more 0 3 is found in the LTP BL than at ACR 
due to pollution from Panama City and occasionally from South American biomass fires. 
(2) The TTL extends lower over LTP than over ACR. It is customary to assume that when 
the free tropospheric segment of an ozone profile resembles BL concentrations that 
convective redistribution has taken place [ Kley et al. 1996; Thompson et al., 1997; Newell 
et al., 1999, Folkins et al., 2006]. Given typical BL concentrations of 20-25 ppbv at the 
two sites, Figure 2 suggests that there are fewer episodes of high-altitude convection over 
LTP than ACR (blue shade, 12-14 km), although that inference may be an artifact of LTP 
having more samples. Consequently, there is less 0 3 in the ACR mean 0 3 profile between 
11-14 km (Figure 3a) than LTP (Figure 3b; cf Figure 13 in Avery et al., 2010), giving a 
slight “S” shape in the ACR 0 3 profile. The “S” is similar to Pacific SHADOZ profiles 
( Folkins et al., 2000; 2006; cf Kley et al., [1996]) and to IONS-06 soundings (< http:// 
croc.gsfc.nasa. gov/in texb/ionso6.html>) over Mexico City (19N, 99 W) that detected the 
highest 0 3 at 8-9 km and minimum 0 3 at 12-14 km ( Thompson et al., 2008; Fig 5). 

Convective signatures in DC-8 sampling, registered as locally minimum ozone in 
profiles from the FASTOZ instrument, were concentrated between 10 and 11 km [ Avery et 
al., 2010]. Elevated lower tropospheric tracers such as methyl hydrogen peroxide and 
organic bromine-containing constituents (Figures 13-15 in Avery et al., 2010) occur 
coincidentally with the lower ozone. Mean DC-8 0 3 profiles and the soundings (Figure 
3a) are somewhat divergent. This is not surprising given that the sondes are at fixed 
locations whereas the DC-8 deliberately sampled in the vicinity of active convection. 

3-2 Overview of Wave Signatures in the Sondes 

Figure 4a shows that the 0 3 labeled GW is concentrated throughout the TTL and 
lower stratosphere over both LTP and ACR, with frequencies similar to those at the 
Paramaribo and San Cristobal SHADOZ sites (Figure 4b). The latter depicts wave 
frequencies averaged for 1999-2007, which are nearly identical to June-July-August 
statistics. The waves in Figures 4a, b are similar to other SHADOZ tropical [ Thompson 
et al., 2010a, b] and northern hemisphere subtropical locations (Figure 4 in Thompson et 
al., 2008). As in Thompson et al. [2010a, b], the GW signal in the TTL is more frequent at 
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the station that is closer to the equator (Las Tablas and San Cristobal, respectively, in 
Figures 4a and 4b). However, ACR has a higher GW frequency at 3 and 6 km, where 
convective outflow often occurs (Section 3.3). Ozone associated with GW is 15-20% of 
the tropospheric column (Table 2). In half the days with TC4 soundings at both LTP and 
ACR, LTP 0 3 at 3-6 km is 40-50 ppbv vs ~30 ppbv over ACR (cf Figure 3 mean 0 3 
profiles). Thus, column tropospheric 0 3 is -30% higher over LTP than ACR (Table 2). 
3.3 Case Studies with Varying Convective and Wave Influence 

Days with coincident launches at LTP and ACR (Figure 5a) offer an opportunity to 
compare ozone profiles at the two sites and to use TC4 aircraft data, satellite imagery and 
meteorological products to interpret convective and wave influences. On all dates 
illustrated in Figure 5a the profiles display convective signatures typically with >10 DU 
of column 0 3 to 20 km designated as GW. The TC4 period started with a phase 
characterized by active convection over the Costa Rica-Panama region and the adjacent 
Pacific (Section 3.3.1); this corresponds to the left-most section relative to the vertical 
dashed line in Figures 5b, c. This is corroborated by satellite and aircraft imagery in 
flights through 22 July 2007 [ Toon et al., 2010]. In a second phase, from 23 July through 

2 August (Section 3.3.2), between the two dashed lines in Figures 5b, c, TC4 aircraft 
continued to sample convection but it was less active at ACR and LTP. The fraction of 
total 0 3 designated GW declined in most soundings (Figures 5b, c), as the GW-affected 
segments in 0 3 profiles retreated from the troposphere to above 17 km. The amount of 
column 0 3 affected by RW increased compared to the pre-24 July period. More active 
convection resumed in a third phase (Section 3.3.3) after 2 August, when targeted 
sampling by the ER-2, DC-8 and WB-57 took place [Toon et al., 2010]. In the following 
sections episodes illustrative of the three phases are examined to verify the link between 
the GW designation and convective activity and the relationship of RW to extra-tropical 
and/or pollution influences. 

3.3.1 13-22 July 2007: Active Convection and Elevated GW 

13 July . Convective activity at LTP and ACR, as given by GW amount to 20 km, is 
similar (Figure 5a). The fraction of FT 0 3 influenced by GW is similar in both cases 
(~20%; Figures 5b, c). A typical pair of profiles (Figure 6a) shows that the ACR 0 3 
concentration averages -40 ppbv from the surface to the tropopause (14-15 km), whereas 

0 3 in the upper troposphere over LTP exceeds 50 ppbv. BL 0 3 at LTP is less than at ACR. 
The RH profiles have considerable structure, with moister air over LTP than ACR. 
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Evidence for convection over ACR comes from the DC -8 flight from California to San Jose, 
Costa Rica, on 13 July 2007 . The last 100-150 km of the flight encompassed a descent 
near ACR at ~ 2 ioo UTC, after the aircraft had crossed the ITCZ. The DC -8 uv-DIAL 
image of 0 3 (Figure 6b) captures the morphology of convective impact throughout the FT 
and TTL. North of the ITCZ (northern edge at 13 N), the FT was penetrated by pollution 0 3 
(>80 ppbv) and aerosols traced to biomass fires interacting with convection. South of the 
ITCZ, just before descent, FT 0 3 dropped to 40-50 ppbv, except for a localized 0 3 
minimum (< 30 ppbv) around 10 km (Figure 6b), similar to FT 0 3 structure over ACR 
(Figure 6a), and to the DC -8 FASTOZ 0 3 during descent. Note a very thin layer of low 
ozone at 13 km in the uv-DIAL image, corresponding to an altitude with low 0 3 over ACR. 
A similar low-0 3 “bubble” sampled on 17 July 2007 is discussed by Petropavlovskikh et al., 
2010 s. The uv-DIAL aerosols (not shown) corresponding to Figure 6b indicate a “clean” 
FT south of the ITCZ except for thin cirrus at the tropopause, consistent with the 
soundings. Convective indicators, elevated CO, methyl-hydrogen peroxide, lightning NO, 
ultrafine particles (not shown; see flight report at < h 1 1 p : / / e s p o . n a s a . go v / 1 c 4 / d o c s > ) , 
penetrated south of the ITCZ. These pollutants came from north of the ITCZ. 
Interhemispheric transport during convection is well known over the Atlantic [Jonquieres 
et al., 1998 ; Thompson et al., 2000 ; Edwards et al., 2003 ]. 

On 13 July (Figure 5b, c) about half the tropospheric 0 3 over both LTP and ACR 
corresponds to RW. The corresponding segments in the ACR profile (Figure 6a) are 
locally enhanced in 0 3 and reduced in RH, suggesting extra-tropical influence. This is a 
good example of a sounding that captures both advective and convective signals. 

22 July. On this day there is more 0 3 in the TTL and LS identified as GW over 
LTP than ACR (Figures 5a, 7a). There is relatively little tropospheric 0 3 associated with 
GW over either one (Figures 5b, c). A convective contrast between LTP and ACR is 
evident in satellite imagery (see ER -2 flight report for 22 July 07 at < http://espo. 
nasa.gov/tc 4 /d 0 cs>) . When the ER -2 sampled near LTP on 22 July, the cloud and 
precipitation CPL-CRS imagery (Figure 7b) indicated several levels of convective 
outflow. For example, there is cirrus with a 13 . 5-14 km cloud top, coinciding with a GW- 
labeled segment in the LTP sounding, right above a localized 0 3 minimum (black profile, 
Figure 7a). The CPL-CRS indicates another cloud outflow layer at 4 km. This region is 
not designated as GW although there is a local 0 3 minimum at 4 km and relatively high 
RH. From 2-5 km the designation is RW; 0 3 is 60 % greater than at the surface, suggesting 
possible pollution transport. Surface ozone at LTP, measured in NATIVE on 22 July, 
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averaged 15 ppbv (Figure 8a) but the moderately high CO mixing ratio, 135 ppbv, 
supports an interpretation of pollution (Figure 8b). 

3.3.2. 23 July-2 August: Elevated RW 

During the period 23-28 July most sonde launches were at LTP. At both sites, the 
soundings from 23 July to 2 August displayed fractionally higher RW 0 3 segments in the 
troposphere (Figure 5b, c). An example appears in Figure 9a. The RW signal over LTP 
extends from 2 to 17 km on 2 August, except for a 2-km region, and there is no GW 
segment. ACR displays convective influence in terms of GW only from 11-15 km. An RW 
segment from 5-10 km coincides with a dry layer from 4-9 km; this is suggestive of air 
from the extra-tropics. Even though convection picks up on 3 August, remnants of the RW 
feature persist over both sites on that day (Figure 9b). 

During the period of greater RW, there were also changes in surface 0 3 and CO at 
Las Tablas (see NATIVE data in Figure 8). On 23 July a normal diurnal 0 3 cycle was 
observed, with a near-zero nocturnal minimum. LID analysis for the 23 July sounding at 
LTP was not valid, indicating active transition and no stable layers. However, on 24 July 
the 0 3 minimum rises to -15 ppbv. This causes daily mean 0 3 to increase from 17 ppbv on 
23 July to 25 ppbv on 24 July. At this time, CO dropped below 90 ppbv (Figure 8b), one 
of the lowest values at NATIVE during TC4. 

3. .3. .3 3-iQ August: Return of Active Convection. Elevated GW 

a August. The fractions of GW and RW tropospheric segments over Las Tablas 
and Alajuela are nearly the same (Figures 5b, c) and the wave structures are similar 
(Figure 9b). Over LTP the GW signal that was absent on 2 August returns at 13-20 km 
(Figure 9b). The DC-8 and ER-2 sampled not far from LTP near active convection (refer 
to GOES image with flight tracks for all three aircraft on 3 August; Figure 16 in Toon et al. 
[2010]). The DC-8 profiles in Figure 10a were taken over NATIVE from 1705-1735 UTC 
during the convective period that was sampled less than 50 km away by the ER-2; the LTP 
launch took place at 1741 UTC. The upper of two 0 3 minima over LTP (10-15 km, Figure 
9b), corresponds to cloud outflow recorded by satellite imagery (not shown) and to cirrus 
in the ER-2 CPL-CRS product (Figure 10b) at 13-15 km. A second 0 3 minimum with 
elevated CO at 5-6 km (Figure 10a) is also located at a cloud outflow layer (lower bar 
arrow in Figure 10b). Over ACR, 0 3 mixing ratios averaged > 80 ppbv in the 6-13 km 
segment (Figure 9b). Above 8 km, the RH drops off sharply, suggesting extra-tropical 
air. The elevated 0 3 coincides with a mostly RW segment that brackets the tropopause, 
which is 1.5 km higher over LTP than ACR (Figure 2). 
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5 August . During the period 3-5 August there was a sharp transition in the LTP 
and ACR profiles that is reflected in the 0 3 budgets (Figures 5a-c); the GW 0 3 budgets to 
20 km are nearly identical over the two stations on 4 August (Figure 5a). The upper 
tropospheric 0 3 mixing ratios decline over both sites, from a mean of -75 ppbv on 3 
August to 45 ppbv on 5 August (Figures 9b, c); during the transition, LTP free 
tropospheric 0 3 declines to 22 DU on 4 August, the lowest 0 3 column in the TC 4 period 
(Figure 5b). From 4 to 5 August there are further transitions in vertical 0 3 structures 
over ACR and LTP (Figure 9c). The GW signal that starts over LTP (ACR) at 12 km (11 
km) on 4 August (not shown) retreats to 15 km on 5 August. Over Las Tablas a several-km 
thick layer that is 20 ppbv above background coincides with RW (Figure 9c; also Figure 
1). DC -8 aircraft sampling near LTP on 5 August confirmed the possibility of extra- 
tropical origins at 8-10 km, where 0 3 mixing ratios were 50-70 ppbv in a relatively dry 
layer [Avery et al., 2010 ]. However, just above this layer, at 11 km, the DC -8 noted cleaner 
than usual conditions, signifying convective outflow of marine boundary layer air (upper 
outflow on Figure 10c). Near the surface, the DC -8 detected pollution due to biomass 
fires (Flight notes at < h U p : / / w w w . e s p o . n a s a . go \ 7 1 c 4 / f 1 i g h t D o c s . p h p > ) . 

The 5 August ozonesonde launch at LTP (1505 UTC) was timed for nearby DC -8 
spiral ("1540 UTC) and ER -2 sampling of convective cells to the southwest (Figure 10c; 
cf Figure 19 in Toon et al., 2010 ). Between 2.5 km and 5.1 km, the sonde bounced up and 
down five times within a convective cell while 0 3 concentrations increased by ~io ppbv. 
Morris et al. [ 2010 ] use satellite imagery, lightning data, radar and OMI N0 2 maps, along 
with DC -8 measurements to demonstrate that lightning NO production is responsible for 
much of the increase. The lower outflows in Figure 10c correspond to the region of the 
bouncing sonde. The profile in Figure 9c is based on the final sonde ascent from 2.5 to 
5.1 km; the GW signal is not detected. 

4 5 _^Vaves s intheJi<|uatorialAmericasj = j99922002^_TC^ 5 ii^Context. 

A perspective for interpreting convective influence over LTP and ACR is provided 
by wave frequencies over Paramaribo and San Cristobal (Figure 4b), SHADOZ sites 
operating since 1999 [ Thompson et al., 2003 a, b]. The amplitudes of individual layers and 
wave structure at Paramaribo and San Cristobal resemble those for LTP and ACR (Figure 
1) as do the structure of the GW frequencies at LTP and ACR (Figures 4a, b). Similar GW 
structure appears over equatorial Indian Ocean sites, eg Watukosek, Kuala Lumpur, that 
display the highest annually averaged GW frequency, ~ 6 o% [ Thompson et al., 2010 a]. 

The higher GW frequency at San Cristobal leads to a larger GW fraction of tropospheric 0 3 
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(Figure ll). Although the tropospheric 0 3 column averaged 25 DU in 2007 , compared to 
28 DU for LTP, 25 % of 0 3 is GW-affected at San Cristobal compared to 15 % at LTP (Table 
2). TC 4 was timed for the onset of the sub-tropical convective season and the North 
American monsoon. Figure 4c, a summary of monthly averaged GW over Paramaribo, 
shows that JJA has about half the maximum GW frequency, typically a December 
occurrence. 

An interannual view of ozone budgets and convective influence appears in Figure 
11, where the mean JJA tropospheric 0 3 column (with segments for RW and GW) is 
displayed for Costa Rica (2006 only), Paramaribo and San Cristobal from 1999 - 2007 . At 
San Cristobal, 2006 is a low-0 3 year compared to the six others, possibly due to a 
moderate El Nino [ Arguez , 2007 ]. In the eastern Pacific, El Nino tends to enhance 
convective activity, mixing lower 0 3 air from the BL throughout the troposphere. [Logan 
et ah, 2008 ]. The GW-affected tropospheric ozone amount in 2006 is only slightly lower 
than normal over San Cristobal but the total tropospheric column dropped from a mean 
22-23 DU ( 1999 - 2005 ; Figure 11) to 18 DU so the GW fraction is magnified. At Heredia 
(20 km away) 2006 tropospheric 0 3 is lower than over ACR in 2007 (Figure 11). 

General meteorological conditions at Paramaribo ( 6 N), Panama ( 8 N), and Alajuela 
( 10 N) are similar, with the ITCZ migrating over each. The SHADOZ wave climatology 
[Loucks, 2007 ; Thompson et ah, 2010 ] shows GW frequency diminishing with increasing 
latitude. A GWI (Gravity Wave Index; Figure 12), based on the fraction of the 0 3 column 
to 20 km denoted as GW, provides a quantitative comparison of site-to-site and 
interannual variability. GWI is larger at San Cristobal than Paramaribo until 2004 when 
data gaps at San Cristobal compromise the record. The gaps also preclude conclusive 
linkage of GWI to signals associated with an El Nino. 

5. Summary 

During TC 4 , in July and early August 2007 , ozonesondes and radiosondes were 
launched several times/week at Alajuela, Costa Rica (ioN, 84 W) to characterize convective 
influences in the tropical tropopause layer (TTL). At Las Tablas, Panama ( 7 . 8 N, 80 W), a 
coastal site 300 km southwest of Panama City, 0 3 profiles from daily sondes, surface 0 3 , 
CO and other tracers are analyzed. Laminar identification provides a systematic approach 
to classifying wave signatures in sounding data, giving a statistical perspective on the TC 4 
period and the longer-term SHADOZ sounding record at Paramaribo and San Cristobal. 
The findings are summarized: 
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453 • GW influences, possibly due to semi-permanent Kelvin waves in the TTL and lower 

454 stratosphere (cf Grant et al., 1998; Thompson et al., 2010a) appeared in 50% 

455 (40%) of Las Tablas (Alajuela) sondes. The GW structure and frequency are similar 

456 to those over SHADOZ stations at San Cristobal and Paramaribo. 

457 • On average there is 35-40% more tropospheric column 0 3 at LTP than ACR during 

458 TC4 and 20% more at LTP than at San Cristobal, a remote marine station, 1400 km 

459 southwest of LTP. 

460 • June- July-August 0 3 budgets at Paramaribo and San Cristobal suggest that 2007 

461 was a “typical” year in terms of tropical equatorial 0 3 amount and convective 

462 activity expressed in GW frequency. During 1999-2006, Paramaribo and San 

463 Cristobal display 0 3 column amounts and convective influence that bracket the TC4 

464 ACR and LTP values. 

465 Classification of wave types through LID is validated through case studies in which 

466 aircraft observations support interpretation of convective influences (with the GW 

467 designation) and extra-tropical impact, corresponding to RW. Laminae of low-0 3 surface 

468 air convectively injected into the free troposphere are detected by LID, frequently 

469 interleaved with the richer-0 3 layers; subtle day-to-day variations are captured. The 

470 pattern of convection inferred from LID is consistent with the meteorological evolution of 

471 the campaign [ Toon et al., 2010]. The early part of TC4, to 22 July 2007, was 

472 characterized by persistent GW throughout the free troposphere and TTL. After a less 

473 active period, from 23 July until approximately 2 August, with lower GW signals and 

474 greater RW, GW increased in frequency along with convection. 

475 Sonde and aircraft data established further the convection-GW linkage and 

476 demonstrated the prevalence of extra-tropical laminae interleaved with layers from 

477 convective outflow throughout the equatorial Americas. In terms of TC4 objectives, our 

478 analysis of ozone structure strengthens the case for convection as a dominant mechanism 

479 for water vapor transport and cirrus formation in the TTL. The persistence of higher-0 3 

480 laminae in the troposphere requires further investigation to determine the extent to which 

481 these layers are remnants of extra-tropical filaments or associated with localized 

482 equatorial waves. 
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Crawford, J. A Al-Saadi). Las Tablas measurements with the NATIVE trailer were assisted by A. Pino and 
L. Jordan (University of Panama); A. M. Bryan and D. Lutz (Valparaiso Univ); Z. Chen and J. L. Tharp 
(PSU). Costa Rican launches were made by UNA students K. Cerna, V. H. Beita, D. Gonzalez. Thanks to 
Mission Scientists M. R. Schoeberl and P. A. Newman for flight notes and to K. E. Pickering for 
discussing lightning data. Thanks to EAB, BvdW, AOG (PSU) for analysis. 

APPENDIX. ABBREVATIONS AND ACRONYMS 

ACP= Alajuela, Costa Rica ( 10 N, 84 W) 

BL = Boundary Layer, here determined from radiosondes (Yorks et al. 2009 ) 

CALIPSO = Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation 
CPL = Cloud Physics Lidar (ER -2 instrument) 

CRS = Cloud Radar System (ER -2 instrument) 

DU = Dobson Unit (1 DU = 2.69 x 10 16 cm' 2 ) 

ECC = Electrochemical Concentration Cell (ozonesonde type) 

FT = Free Troposphere 

GOES = Geostationary Operational Environmental Satellites 
GSFC = Goddard Space Flight Center 
GWI = Gravity Wave Index 

INTEX = Intercontinental Transport Experiment (-A, 2004 ; B, 2006 ) 

IONS = INTEX Ozonesonde Network Study < http://croc.gsfc.nasa.gov/intex/ions.htmI; ... 
intexb/ionso6> 

ITCZ = Intertropical Convergence Zone 

JJA = June-July-August 

LTP = Las Tablas, Panama ( 7 . 8 N, 80 W) 

NASA = National Aeronautics and Space Administration 

OMI = Ozone Monitoring Instrument 

P-T-U = Pressure-Temperature-Humidity 

RH = Relative Humidity 

RWI = Rossby Wave Index 

SHADOZ = Southern Hemisphere Additional Ozonesondes 
<http: / /croc.gsfc. nasa.gov/shadoz> 

SOWER = Stratospheric Ozone and Water Vapor in Equatorial Regions 
TC 4 = Tropical Composition, Clouds and Climate Coupling ( 2007 ) 
< http://espo.nasa.gov/tc4> 

TCSP = Tropical Convective Systems and Processes (Costa Rica, 2005 ) 

TTL = Tropical Tropopause Layer (sometimes tropopause transition layer) 

UV-DIAL = Ultraviolet Differential Airborne Lidar [Laser Detection and Ranging] 

Table 1 . Stations for which data are used. Further technical details given in Table A-i in 
Thompson et al. ( 2003 a) and in Thompson et al. ( 2007 b). 

Station Latitude, Instrument Type, Radiosonde Co-Investigator/ 

Longitude Sensing Solution Sponsor 

Las Tablas 7.8N, 80W ENSCI, 0.5% KI, buffered RS-80-15N G. A. Morris, A. M. Thompson 
Heredia/Alajuela, io.oN,84 W ENSCI, 1% KI, RS-80& H. Vomel; J. Valverde Canossa 

buffered Cryogenic Frost-point hygrometer 

San Cristobal 0.92S, 89.6W ENSCI 2% unbuffered RS-80 H. Vomel. INAM HI (National 

to 2006; 1%, reduced Inst, of Hydrology and 

buffer, 2006- Meteorology of Ecuador), M. V. A. 

Reyes [ Johnson et al., 2002; Thompson 
et al., 2007b] 

Paramaribo 5-8N,55W SPC, 1% KI buffered RS-80 to 2005 G. Verver & Met. Service 

RS-92, 2005- Suriname [ Peters et al., 

2003; Fortuin et al., 2006] 
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Table 2. Free tropospheric ozone columns during June-July-August 


2007. ACR mean omits 28 July sounding where data ended below the tropopause. 


Station 

GW O., 

RWO., 

Other O., 

Total 

ACR - DU 

2.9 

8.2 

9-3 

20.4 

ACR - % 

14 

40 

46 

100 

LTP - DU 

3-94 

13-2 

10.8 

28 

LTP - % 

15 

47 

38 

100 

San Cris.- DU 

5-5 

7-5 

12 

25 

San Cris.- % 

24 

28 

48 

100 


References 

Arguez, A., ed. (2007): State of the climate in 2006. Bull. Am. Meteor. Soc, 88, S1-S135. 

Avery, M. A., J. Joiner, C. Twohy, D. McCabe, E. Atlas, D. Blake, P. Bui, J. Counse, J. Dibb, G. 

Diskin, R. Gao, P. Lawson, M. McGill, D. Rogers, G. Sachse, R. Salawitch, E. Scheuer, 

K. Severance, A. M. Thompson, C. Trepte, P. Wennberg, J. Ziemke, (2010) Convective 
distribution of tropospheric ozone and tracers in the central American ITCZ Region: 
Evidence from observations during TC4, J. Geophys. Res., doi: 10.1029/2009 
JD013450, in press. Manuscript available at TC4 website. Contact 

btoon(a) lasp.colorado.edu for password. 

Bertram, T., et al. (2007), Direct measurements of the convective recycling of the upper 
troposphere, Science, 315, 816-820, doi:io. 1126/science. 1134548. 

Bloom, S., et al. (2005), Documentation and validation of the Goddard Earth Observing System 
(GEOS) data assimilation system - Version 4. Technical Report Series on Global 
Modeling and Data Assimilation 104606. 

Browell, E. V., et al. (1996), Ozone and aerosol distributions and air mass characteristics over 
the South Atlantic Basin during the burning season, J. Geophys. Res., 101, 24,043- 
24,068. 

Cooper, O. R., et al. (2006), Large upper tropospheric ozone enhancements above mid-latitude 
North America during summer: In situ evidence from the IONS and MOZAIC ozone 
networks, J. Geophys. Res., ill, D24S05, doi: 10.1029/2006JD007306. 

Deshler, T., et al. (2008), Balloon experiment to test ECC-ozonesondes from different 

manufacturers, and with different cathode solution strengths: Results of the BESOS 
flight, J. Geophys. Res., 113,004307, doi:io.i029/2007JDoo8975. 

Dougherty, K. M. (2008), The effect of ozonopause placement on tropospheric ozone budgets: 
An analysis of ozonesonde profiles from selected IONS-06 sites, MS Thesis, The 
Pennsylvania State University. 

Draxler, R. R., and G. D. Rolph (2003), HYSPLIT (Hybrid Single-Particle Lagrangian 

Integrated Trajectory) model, http://www.arl.n0aa.gov/ready/hysplit4.html, NOAA 
Air Resour. Lab., Silver Spring, MD. 

Edwards, D. P., et al. (2003), Tropospheric ozone over the tropical Atlantic: A satellite 
perspective, J. Geophys. Res., 108, 4237, doi: 10.1029/2002JD002927. 

Fast, J. D., et al. (2007), A meteorological overview of the MILAGRO field campaign, Atmos. 
Chem. Phys. 7, 2233-2257. 

Folkins, I., S. J. Oltmans, and A. M. Thompson (2000), Tropical convective outflow and near- 
surface equivalent potential temperatures, Geophys. Res. Lett., 27, 2549-2552. 


17 


584 

585 

586 

587 

588 

589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 

600 

601 

602 

603 

604 

605 

606 

607 

608 

609 

610 

611 

612 

613 

614 

615 

616 

617 

618 

619 

620 

621 

622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

632 

633 

634 

635 


Folkins, I., P. Bernath, C. Boone, K. Walker, A. M. Thompson, and J. C. Witte (2006), The 

seasonal cycles of O3, CO and convective outflow at the tropical tropopause, Geophys. 
Res. Lett., 33, L16802, doi:io.i029/2C>o6GLo266o2. 

Fortuin, P., et al. (2007), Origin and transport of tropical cirrus clouds observed over 
Paramaribo, Suriname (5.8°N, 55.2°W), J. Geophys. Res., 112 , D09107, 
doi:io. 1029/2005 JD006420. 

Fueglistaler, S., A. E. Dessler, T. J. Dunkerton, I. Folkins, Q. Fu, and P. W. Mote (2009), 

Tropical tropopause layer, Rev. Geophys., 47 , RG1004, doi:io.i029/20o8RGooo267. 

Fujiwara, M., K. Kita, and T. Ogawa (1998), Stratosphere-troposphere exchange of ozone 
associated with the equatorial Kelvin wave as observed with ozonesondes and 
rawinsondes, J. Geophys. Res., 103 , No. D15, 19,173-19,182. 

Fujiwara, M., F. Ilasebe, M. Shiotani, N. Nishi, II. Vomel, and S. J. Oltmans (2001), Water 
vapor control at the tropopause by equatorial Kelvin waves observed over the 
Galapagos, J. Geophys. Res., 28 , 3143-3146. 

Grant, W. B., R. B. Pierce, S. J. Oltmans, and E. V. Browell (1998), Seasonal evolution of total 
and gravity wave induced laminae in ozonesonde data in the tropics and subtropics, 
Geophys. Res. Lett. 25, 1863-1866. 

Hasebe, F., M. Fujiwara, N. Nishi, M. Shiotani, II. Vomel, S. Oltmans, II. Takashima, S. 

Saraspriya, N. Komala, and Y. Inai (2007), In situ observations of dehydrated air 
parcels advected horizontally in the Tropical Tropopause Layer of the western Pacific, 
Atmos. Chem. Phys., 7, 803-813. 

Hlavka, D., L. Tian, W. Hart, L. Li, M. McGill, and G. Ileymsfield (2010), Vertical cloud 

climatology during TC4 derived from high-altitude aircraft merged lidar and radar, J. 
Geophys. Res., this issue. 

Johnson, B. J., S. J. Oltmans, H. Vomel, T. Deshler, C. Kroger, and H. G. J. Smit (2002), ECC 
ozonesondes pump efficiency measurements and sensitivity tests of buffered and 
unbuffered sensor solutions, J. Geophys. Res., ic> 7 ( D19), 4393, doi: 
10.1029/2001JD000557. 

Jonquieres, I., A. Marenco, A. Maalej, and F. Rohrer (1998), Study of ozone formation and 
transatlantic transport from biomass burning emissions over West Africa during the 
airborne Tropospheric Ozone Campaigns TROPOZ I and TROPOZ II, J. Geophys. Res., 
103 , 19,059-19,073. 

Kley, D., P. J. Crutzen, H. G. J. Smit, H. Vomel, S. J. Oltmans, H. Grassl and V. Ramanathan 
(1996) Observations of near-zero ozone over the convective Pacific: Effects on air 
chemistry, Science, 274 , 230-233. 

Logan, J. A., I. Megretskaia, R. Nassar, L.T. Murray, L. Zhang, K.W. Bowman, II. M. Worden, 
and M. Luo (2008), Effects of the 2006 El Nino on tropospheric composition as 
revealed by data from the TES, Geophys. Res. Lett., 35 , L03816, 
doi:io.i029/2007GLo3i698. 

Loucks, A. L. (2007), Evaluation of dynamical sources of ozone laminae in the tropical 
troposphere and tropical tropopause layer, M.S. Thesis, Penn State University. 

McGill, M. J., L. Li, W. D. Hart, G. M. Heymsfield, D. L. Hlavka, P. E. Racette, L. Tian, M. A. 
Vaughan, and D. M. Winker (2004), Combined lidar-radar remote sensing: Initial 
results from CRYSTAL-FACE, J. Geophys. Res., 109 , D07203, 
doi:io.i029/2003JDc>04030. 

Morris, G. A., et al. (2010), Observations of ozone production in a dissipating convective cell 
during TC4, J. Geophys. Res., doi: 10.1029/2009JD 014098, submitted. Manuscript 
available at TC4 website. Contact btoon@lasp.colorado.edu for password. 

Newell, R. N., V. Thouret, J. Y. N. Cho, P. Stoller, A. Marenco, and H. G. Smit (1999), Ubiquity 
of quasi-horizontal layers in the troposphere, Nature, 398 , 316-319. 

Oltmans, S.J., et al. (2001), Ozone in the Pacific tropical troposphere from ozonesonde 
observations, J. Geophys. Res., 106 , 32503-32526. 


18 


636 

637 

638 

639 

640 

641 

642 

643 

644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

658 

659 

660 

661 

662 

663 

664 

665 

666 

667 

668 

669 

670 

671 

672 

673 

674 

675 

676 

677 

678 

679 

680 

681 

682 

683 

684 

685 

686 

687 

688 


Oltmans, S. J., et al. (2004), Tropospheric ozone over the North Pacific from ozonesonde 
observations, J. Geophys. Res., 109, D15S01, doi: 10.1029/2003JD003466. 

Peters, W., P. Fortuin, H. Kelder, C.R. Becker, J. Lelieveld, P.J. Crutzen, and A.M. Thompson 
(2004), Tropospheric ozone over a tropical Atlantic station in the Northern 
Hemisphere: Paramaribo, Surinam (6°N, 55°W), Tellus, 56, 21-34. 

Petropavloskikh, I., E. Ray, S. M. Davis, K. Rosenlof, G. Manney, R. Shetter, S. Hall, K. 

Ullmann, L. Pfister, J. Hair, M. Fenn, M. Avery, and A. M. Thompson, (2010) Low ozone 
bubbles observed in the tropical tropopause layer during the TC4 campaign in 2007, J. 
Geophys. Res., doi: 10.1029/2009JD012804, in press. Manuscript available at TC4 
website. Contact btoon(S lasp.colorado.edu for password. 

Pierce, R. B., and W.B. Grant (1998), Seasonal evolution of Rossby and gravity wave induced 
laminae in ozonesonde data obtained from Wallops Island, Virginia, Geophys. Res. 

Lett., 25, 1859-1862. 

Piotrowicz, S. R., H. Bezdek, G. Harvey, and M. Springer-Young (1991), On the ozone minimum 
over the equatorial Pacific Ocean. J. Geophys. Res. 96, 18679-18687. 

Randel, W. J., D. J. Seidel, and L. L. Pan (2007), Observational characteristics of double 
tropopauses, J. Geophys. Res., 112, D07309, doi:io.i029/20o6JDoo7904. 

Randriambelo, T., J-L. Baray, S. Baldy, A. M. Thompson, S. J. Oltmans, and P. Keckhut (2003), 
Investigation of the short-term variability of tropical tropospheric ozone, Annales 
Geophysiques, 21, 2095-2106. 

Read, K. A. et al. (2008), Extensive halogen-mediated ozone destruction over the tropical 
Atlantic Ocean, Nature, 453, 1232-1235, doi: I0.i038/nature07035. 

Schoeberl, M. R., and L.C. Sparling (1995), Trajectory modeling: Diagnostic tools in 

atmospheric physics, S. I. F. Course CXVI, edited by G. Fiocco and C. Visconti, North- 
Holland, Amsterdam. 

Schoeberl, M. R., et al. (2006), Overview of the EOS aura mission, IEEE Trans., 44 (5), 1066- 
1074, doi:io.no9/TGRS. 2005. 861950. 

Selkirk, H. B., H. Vomel, J. M. Valverde Canossa, L. Pfister, J. A. Diaz, W. Fernandez, J. 

Amador, W. Stolz, and G. Peng, The detailed structure of the tropical upper troposphere 
and lower stratosphere as revealed by balloonsonde observations of water vapor, ozone, 
temperature and winds during the NASA TCSP and TC4 Campaigns, J. Geophys. Res., 
in press. ** Manuscript available at TC4 website. Contact btoon( 5 >lasp.colorado.edu for 
password. 

Smit, H. G. J., et al. (2007), Assessment of the performance of ECC-ozonesondes under quasi- 
flight conditions in the environmental simulation chamber: Insights from the Jiilich 
Ozone Sonde Intercomparison Experiment (JOSIE), J. Geophys. Res., 112, D19306, doi: 
10. 1029/2006 JD007308. 

Takashima, H., and M. Shiotani (2007), Ozone variation in the tropical tropopause layer as 
seen from ozonesonde data, J. Geophys. Res., 112, D11123, 
doi:io.i029/20o6JDoo8322. 

Teitelbaum, H., J. Ovarlez, H. Kelder, and F. Lott (1994), Some observations of gravity-wave- 
induced structure in ozone and water vapour during EASOE, Geophys. Res. Lett., 21, 
1483-1486. 

Teitelbaum, H., et al. (1996), The role of atmospheric waves in the laminated structure of ozone 
profiles at high latitude. Tellus, 48A, 442-455. 

Thompson, A. M., et al. (1993), SAGA-3 ozone observations and a photochemical model 

analysis of the marine boundary layer during SAGA-3, J. Geophys. Res., 98, 16955- 
16968. 

Thompson, A. M., et al. (1996), Where did tropospheric ozone over southern Africa and the 

tropical Atlantic come from in October 1992? Insights from TOMS, GTE/TRACE-A and 
SAFARI-92, J. Geophys. Res., 101, 24,251-24,278. 

Thompson, A. M., W.-K. Tao, K. E. Pickering, J. R. Scala, and J. Simpson (1997), Tropical deep 
convection and ozone formation, Bull. Amer. Met. Soc., 78, 1,043-1,054. 


19 


689 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 

701 

702 

703 

704 

705 

706 

707 

708 

709 

710 

711 

712 

713 

714 

715 

716 

717 

718 

719 

720 

721 

722 

723 

724 

725 

726 

727 

728 


Thompson, A. M., B. G. Doddridge, J. C. Witte, R. D. Hudson, W. T. Luke, J. E. Johnson, B. J. 
Johnson, S. J. Oltmans, and R. Weller (2000), A tropical Atlantic paradox: Shipboard 
and satellite views of a tropospheric ozone maximum and wave-one in January- 
February 1999, Geophys. Res. Lett., 27,3317-3320. 

Thompson, A. M., et al. (2003a), Southern Hemisphere Additional Ozonesondes (SIIADOZ) 
1998-2000 tropical ozone climatology. 1. Comparison with TOMS and ground-based 
measurements, J. Geophys. Res., 108 , 8238, doi: 10.1029/2001JD000967. 

Thompson, A. M., et al. (2003b), Southern Hemisphere Additional Ozonesondes (SIIADOZ) 
1998-2000 tropical ozone climatology. 2. Tropospheric variability and the zonal wave- 
one, J. Geophys. Res., 108 , 8241, doi: 10.1029/2002JD002241. 

Thompson, A. M., et al. (2007a), IONS (INTEX Ozonesonde Network Study, 2004). 1. 

Summertime UT/LS (Upper Troposphere/Lower Stratosphere) ozone over northeastern 
North America, J. Geophys. Res., 112 , D12S12, doi: 10.1029/2006JD007441. 

Thompson, A. M., J. C. Witte, II. G. J. Smit, S. J. Oltmans, B. J. Johnson, V. W. J. II. Kirchhoff, 
and F. J. Schmidlin (2007b), Southern Hemisphere Additional Ozonesondes (SIIADOZ) 
1998-2004 tropical ozone climatology. 3. Instrumentation, station variability, 
evaluation with simulated flight profiles, J. Geophys. Res., 112 , D03304, doi: 10.1029/ 
2005JD007042. 

Thompson, A. M., J. E. Yorks, S. K. Miller, J. C. Witte, K. M. Dougherty, G. A. Morris, D. 

Baumgardner, L. Ladino, and B. Rappenglueck (2008), Tropospheric ozone sources and 
wave activity over Mexico City and Houston during Milagro/Intercontinental Transport 
Experiment (INTEX-B) Ozonesonde Network Study, 2006 (IONS-06), Atmos. Chem. 

Pf jys., 8,5113-5125. 

Thompson, A. M., A. L. Loucks, S. Lee, S. K. Miller (2010a), Gravity and Rossby wave 

influences in the tropical troposphere and lower stratosphere based on SHADOZ 
(Southern Hemisphere Additional Ozonesondes) soundings, 1998-2007, J. Geophys. 
Res., doi: 10 . 1029 / 2009 JD 013429 , submitted. 

Thompson, A. M., S. J. Oltmans, D. W. Tarasick, P. Von der Gathen, II. G. J. Smit, J. C. Witte 
(2010b), Strategic ozone sounding networks: Review of design and accomplishments, 
Atmos. Environ., in press. 

Toon, O. B., et al. (2010), Planning and implementation of the Tropical Composition, Cloud and 
Climate Coupling Experiment (TC4), J. Geophys. Res., this issue. Manuscript available 
at TC 4 website. Contact btoon@lasp.colorado.edu for password. 

Torres, A. L., and A. M. Thompson (1993), Nitric oxide in the equatorial Pacific boundary layer: 
SAGA-3 measurements, J. Geophys. Res., 98 , 16949-16954. 

Yorks, J. E., A. M. Thompson, E. Joseph, and S. K. Miller (2009), The variability of free 

tropospheric ozone over Beltsville, Maryland (39N, 77W) in the summers 2004-2007, 
Atmos. Environ., 43 , 1827-1838. 

Zafiriou, O. C., M. McFarland, and R. H. Bromund (1980) Nitric oxide in seawater, 

Science, 207 , 637-639, 1980. 


FIGURE CAPTIONS - 

Fig 1 Application of laminar identification (LID) method to typical sounding from Panama. 

Illustrated are normalized 0 3 (solid line), potential temperature (dotted line) and 
correlation between the two quantities (dashed). Correlation criteria for Rossby waves 
(RW) are within vertical lines between - 0.3 and + 0.3 (light blue). The latter designation is 
used in discussion of profiles and budgets. Gravity wave (GW) criterion of Pierce and 
Grant ( 1998 ; see their Figure 1 ) and Thompson et al. ( 2007 a; Figure 3 ) calls for 
normalized 0 3 and 0 correlation to reach 0.7 (vertical line; light green for budgets). For 
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computation of the GW Index, a more restrictive criterion is used, namely, the 
corresponding 0 3 layer amplitude must exceed o.i (10%), as in the darker green. An RW 
Index, not used here, is based on counting only ozone within dark blue. 

Fig 2 Curtain plots of 0 3 mixing ratio to 18 km during TC4 over (a) Alajuela, Costa Rica (ACR); 

(b) Las Tablas, Panama (LTP). White dots refer to the ozonopause as described in the text. 

Fig 3 Mean profiles of ozone, temperature, relative humidity (RH) from surface to 20 km for: (a) 

Alajuela, Costa Rica (ACR); eight 0 3 profiles with slight interference from volcanic S 0 2 
have been smoothed at 3 km. For comparison, mean DC-8 profiles from the FASTOZ 
ozone instrument (green dots) are displayed. The profiles include landing and takeoffs 
from San Jose airport near Alajuela. (b) Las Tablas, Panama. In the latter case, seven 
questionable RH profile segments are omitted from mean. 

Fig 4 (a) Frequency of GW occurrence over LTP, ACR during July-August 2007 TC4 sampling; 

(b) mean GW frequency over Paramaribo and San Cristobal, based on all 1999-2007 
profiles. For (b) J-J-A and full year means are virtually the same because J-J-A 
frequencies fall about halfway between the annual maximum and minimum frequencies. 
Paramaribo did not launch during TC4; (c) annual GW frequency at Paramaribo. The 
latter is typical of near-equatorial SHADOZ sites [ Thompson et ah, 2010a]. 

Fig 5 (a) Amounts of 0 3 (in DU) from top of the BL to 20 km, affected by GW, RW determined 

by LID [ Thompson et al., 2007a] based on 0 3 and P-T-U soundings from days with both 
Las Tablas (LTP) and Alajuela (ACR) launches during TC4. (b) Same as (a) except for free 
tropospheric 0 3 segment of all LTP soundings during TC4. The free troposphere is defined 
from the top of the BL to the ozonopause, as illustrated in Figure 2; (c) same as (b) for 0 3 
over ACR. The vertical dashed line distinguishes phases of convective activity (greater 
before 23 July 2007 and after 2 August, diminished in between), as detected in the sondes. 

Fig 6 (a) Ozone, RH, temperature profiles at ACR and LTP for 13 July 2007. Vertical bars refer to RW 

(blue) and GW (green) as described in Figure 1. (b) Uv-DIAL image of ozone from DC-8 
flight from California to Costa Rica. Ozone < 40 ppbv, purple, is near surface and also at 
cloud outflow level, ~ 10 km, south of the ITCZ; the latter is the cloudy region at 1945 UTC. 
Note a thin cloud-outflow ozone lamina at 13 km. 

Fig 7 (a) Ozone, RH, temperature profiles at ACR and LTP for 22 July 2007. (b) convective cells with 
outflow layers denoted by arrows on 22 July 2007 near LTP from ER-2 Cloud Physics 
Lidar and Cloud Radar System composite image [ Hlavka et ah, 2010]. 

Fig 8 Daily mean mixing ratios of (a) 0 3 ; (b) CO, from NATIVE in Las Tablas, Panama (7.8N,8o.oN). 

NATIVE CO readings are higher than the median DC-8 profile data (Figure 2 in Avery et 
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ah, 2010), although NATIVE NO/NO y and S 0 2 suggest that pollution is infrequent, 

Fig 9 Profiles from August 2007 TC4 sampling (a) 2 August 2007; (b) 3 August; (c) 5 August. 

Labels as in Figure 6a. For 5 August, the DC-8 0 3 measurement from profiling near LTP 
displayed a high- 0 3 , low-CO layer [Avery et ah, 2010] at 8-10 km. Also on 5 August, the 
ozonesonde package, caught in dissipating convection, was buffeted in updrafts and 
downdrafts between 2 and 5 km, presumably due to balloon icing [Morris et al., 2010]. 
Only the final ascent profile appears in (c), so a GW signal indicating convection does not 
appear below 7 km. However, the RW segment may denote ozone from pollution. 

Fig 10 (a) Ozone, CO from DC-8 spiral on 3 August 2007 at 1505-1535 UTC, suggests extra- 

tropical influence at 6-8 km, with convection at 5 km and above 8 km. (b) 3 Aug 2007 ER- 
2 sampling produced composite CPL-CRS image with convective cells (outflow at 
horizontal arrows) at 7.5N, 80. 5W. Vertical arrows mark LTP sonde launch, (c) Same as 
(b) except for 5 August. 

Fig 11 Averaged ozone amounts (in DU) in the free troposphere affected by GW, RW determined 

by the laminar method using all 0 3 and P-T-U sounding from J-J-A 1999-2007 for San 
Cristobal, Galapagos (o.89S,9oW) and Paramaribo, Suriname (5.8N,55E) and, since 2005, 
for two Costa Rican launch sites near San Jose (Heredia in 2005-2006; Alajuela, ~20 km 
distant, in 2007). The 2007 data at Las Tablas, Panama, are from TC4. BL 0 3 amounts 
(not shown) are 2 DU at San Cristobal, 3.5 DU at Paramaribo and the Costa Rican sites. 

Fig 12 Gravity wave Indices (GWI) based on 0 3 and P-T-U soundings from the Paramaribo and San 
Cristobal SHADOZ sites. 
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